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FOREWORD 


The Industrial Mathematics Society is a professional organiza- 
tion whose objective it is to extend the understanding and application 
of mathematics in industry. Founded early in 1949, it marked the 
first organized attempt of any group to foster a closer relationship 
between mathematics and industry. 


It promotes the cause of mathematics in industry through 
monthly meetings at which technical papers are presented, through 
panel groups which concentrate in specialized fields of knowledge, 
and through publication of worthwhile papers in its annual volume. 


Membership is open to engineers, scientists, and mathematicians 
who are concerned with the effective use of mathematics in industry. 
Annual dues are $5.00 per year including a subscription to the an- 
nual volume, Industrial Mathematics. Application forms and infor- 
mation may be obtained by writing to the Secretary, Industrial 
Mathematics Society, 100 Farnsworth, Detroit 2, Michigan. 
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The Nature of the Helical Gear Oblique 
Contact Line 


By E. J. Wellaver 
Director, Research and Development 
The Falk Corporation 


The oblique line-of-contact of helical gears has certain charac- 
teristics which merit attention by the practical gear designer as 
well as by the mathematical analyst. 

The “oblique” line-of-contact of helical gears in contrast to the 
*“Parallel-to-axis” contact of spur gears provides a third degree of 
freedom. An understanding of this greater versatility in the nature 
of the helical gear oblique contact line leads to its useful control for 
optimum gear performance. 

This paper abstracts a lecture presented before the Industrial 
Mathematics Society. 

As a means to improve the smoothness of spur gear action the 
face width is sometimes divided and the individual sections stag- 
gered to develop step gears. With continued divisions, the appear - 
ance of the helical gear tooth is approached as shown in Figure 1. 
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It is important to indicate that the ultimate contact does not become 
merely an extension of spur gear action, but develops many special- 
ized characteristics. 

The nature of the oblique contact line shall be reviewed as it 
affects the following: 


Contact action 

Total length of oblique contact lines 
Stress pattern 

Load distribution 

Elastic characteristics 
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1. The Contact Action 


There is a characteristic difference between the mating of spur 
and helical teeth particularly in the transmission and transfer of 
load as the contact progresses along the line of action. 
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ZONE OF ACTION -TWO TOOTH CONTACT 
Fig. 2. Contact of Spur Gears. 


Tooth action of spur gears normally results in either one or two 
teeth in contact. There can be no axial overlap of the lines of con- 
tact parallel to the axis. These characteristics are illustrated in 
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Fig. 3. Contact of Helical Gears. 
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HELICAL GEAR OBLIQUE CONTACT LINE 3 


Figure 2. In contrast, the helical gear usually has a multiple num- 
ber of teeth in contact depending upon the geometric proportions. 
The lines-of-contact are oblique to the axis as shown in Figure 3. 
The geometric location of the oblique line is shown in Figure 4. 


a 





Fig. 4. Inclination Angle, V of Oblique Contact Line. 


Spur Tooth Load Transmission 





The theoretical transmission and transfer of load for a spur 
gear is characterized by the load being shared alternately by one or 
two teeth for discrete portions of the tooth profile. Without proper 
accuracy load sharing between teeth is not possible and the load is 
abruptly transferred from one tooth to the next. Only one tooth is 
then in contact at any instance of action. Table 1 shows the accuracy 
required to insure load sharing. 

The elastic flexibility is continuously changing throughout the 
contact cycle both because of the location of the load on the profile 
and also due to the fact that either one or two teeth are in contact. 

Any error of spacing directly affects the uniformity of the angu- 
lar velocity transmitted. 


Helical Tooth Transmission 





With a helical gear, profile contact can occur simultaneously 
along an oblique line extending from the beginning to the end of con- 
tact along the line of action. 
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Table I 


Limiting Error in Action for Steel Spur Gears 
(Variation in Base Pitch) 

















| Number | Allowable Error When Amount of Error When 
of | Teeth Share Load Teeth Fail to Share Load 
aay L Load Per In. of Face Load Per In. of Face ] 
sy | 500 lb. | 1,000 lb. | 2,000 lb. 500 lb. | 1,000 lb. | 2,000 Ib.| 
| 15 | 0.0004 0.0007 0.0014 0.0006 0.0011 0.0023 | 
| 20 | 0.0003 | 0.0006 | 0.0011 | 0.0006 | 0.0011 | 0.0023 | 
| 25 | 0.0002 | 0.0005 0.0009 0.0006 | 0.0011 0.0023 

















NOTE: When errors lie between the values indicated, more detailed studies or 
tests are required to determine the load sharing characteristics. 


The load can be made to be shared between a sufficient number 
of teeth to allow the load to be more easily picked up by the incom- 
ing tooth and dropped by the outgoing tooth. 

The elastic flexibility does not change measurably during the 
cycle of contact. 

Spacing errors are less likely to affect uniform angular velocity 
because many phases of profile contact are simultaneously oper- 
ating. 

Conditions of abrasive wear do not cause excessive changes in 
the conjugate profile because accelerated wear along any portion of 
the oblique line-of-contact merely transfers the loading to another 
portion. Spur gears tend to depart from their conjugate profile 
more and more as abrasive wear or lapping proceeds because the 
highest lapping tendencies are concentrated in the position of the 
highest sliding. It is for this reason that helical gears can be suc- 
cessfully lapped and the statement is made “that helical gears can 
‘wear in’ before ‘wearing out’.” 

The total length of the contact developed between mating helical 
teeth is subjected to the effect of more variables than that of spur 
gears. A knowledge of these variables enables advantageous effects 
to be secured. 


2. Total Length of Oblique Contact Lines 


For accurate spur teeth, the total length of contact between mat- 
ing profiles equals the face width when one tooth is in contact and 
twice the face width when two teeth are in contact. 

For helical teeth the total length of contact along the oblique 
lines is a complex relationship between geometric dimensions, 
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particularly the ratio of face width to axial pitch (face contact ratio, 
my ) and the ratio of the length of action in the transverse plane to 
the base pitch (profile contact ratio, m, )- 

The minimum length of contact occurs when one of the oblique 
lines intersects and is external to the contact zone at the outer edge 
as shown at points W and Z in Fig. 5a. The maximum contact oc- 
curs when one of the contact lines intersects from the interior as 
shown at points K and Y, Fig. 5a. 

The result is a discontinuous relationship which can be most 
easily solved by a criteria determination of certain parameters de- 
pending upon the fractional portions of the face contact ratio mp and 
the profile contact ratio m, . 


Total Length of the Oblique Lines of Contact 





NOTE: Attention is directed to these special notations which will 
be used in the following formulas: 


K, and n, are the whole number and fractional portion re- 
spectively of the face contact ratio (axial overlap). For 


example, mr = “% = 4.85 then K, = 4.0 and n, = 0.85. 

x 
K, and n, are the whole number and fractional portion re- 
spectively of the profile contact ratio. For example, mp 


~~ = 1.32 then K, = 1.0 and n,; = 0.32. 


b 
K and K,,, are the number of contact lines intersecting side 
K-Z of the contact zone (Fig. 5a) when a contact line inter- 
sects W or Y for L,,;, and L,,,,, respectively. 


Limin Calculation 


The minimum total length of contact is obviously important be- 
cause this is the condition which develops the maximum load in- 
tensity along the oblique contact lines and governs the profile dura- 
bility (pitting resistance) and sometimes also the maximum bending 
stress. Its ratio to the maximum total length of contact has an im- 
portant bearing on load capacity, mesh stiffness, etc. 


The general equation for Lyin is: 


_ (K, - K, + K)(Z - Kp,) _ KF 


oe sin), cos Vy) 





Criteria equations useful for determining the value of K can be 
derived from Fig. 5b. The change from case I (K = K.) to case II 
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Fig. 5. 


(K = K+ 1) occurs when the angle a o b equals and becomes greater 
than angle aoc. The change takes place when 


(1 - n.) p/n, p, = tany, 





or 
1 - Rn, 
=] 
n, 
Hence for: 
Case I 
L « B. 
-21 then K = K, 
Na 
and 
= (ZF Py ) = n,m, Py 
min 2 
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Jase II 
to > 
- 1 then K = K, + 1 
n 
a 
and 
(ZF/p,) - (1 - n,)(1 - n,) p, 
Li, in — cos y 3 
b 
L Calculation 


max 





The maximum total length of the lines-of-contact is important 
because this often is the condition when the position of the worst 
tooth loading is secured in regard to bending strength and also is 
useful as a measure of the variation existing in the contact cycle. 
The general equation for L,,,,, is: 

_ (Ka + K, + 1 - Ky, )(Z - Kp Pp) " Kin F 


max ~ 


L 





siny, cos VY; 


Criteria equations for the determination of K,,, can be derived from 
Fig. 5c similar to the procedure used for the criteria equations 
for K. 
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Case A 
n. =n, then K K . 
and 
L (ZF/p, ) + Ry (1 - ng) Px 5 
max — cos 
b 
Case B 
n, n, then K,, = K 1 
and 
(ZF/p,) + n, (1 - n_) 
L max = Po S r- Px 6 





cosy, 


The mathematical relations for L minand L max existing through 
one base pitch movement of tooth contact along the length of the line 
of action in the transverse plane is illustrated by Fig. 6 and sum- 
marized for all cases in Fig. 7. 


When either n, or n, is zero, in other words, if the ratios mr 


a BASE PITCH- Py = 1.0 ——» 


CASE Af=-{n,- n> 
CASE B (Nn, -Ng =} 
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Fig. 7. Duration of L Through One Base Pitch Cycle. 
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[case e| [case a] 
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jm AXIAL PITCH ( Px) — eee re 


L 
Fig. 8. “F Ratio vs Axial Overlap. 


or mp, are whole numbers, then the contact length remains un- 
changed during all phases of the tooth contact cycle and is equal to 


L=2Z 
Px 


If the total contact length is integrated in relation to the period of 
duration (Fig. 7) then L,,, is equal to ZF/p,. By general usage 
this has become known as the average contact length rather than the 
arithmetical average of the maximum and minimum contact lengths. 

The ratio L/F as the face increases from one integral axial 
overlap to the next highest is illustrated in Fig. 8. The effect of in- 
creasing face contact ratio (axial overlap) is illustrated in Fig. 9. 
Of significance in Fig. 9 is the fact that the ratio L,,,,/Lmin can 
be high when F/p, is less than one but becomes progressively lower 
with increasing numbers of axial overlaps. The desirability of 
maintaining face contact ratios near an integer to minimize total 
contact length variations is evident. 

This information has practical application in design as can be 
shown by the proportions compared in Table II. It is obvious that 
the changes in Lyin are relatively small and the influence of pitch 
or number of teeth in the pinion is minor for the low ratio gear set 
of Example 1. When low numbers of teeth in the pinion are secured 
with the coarse pitches shown in Example 2, an appreciable loss 
in contact can be present. For the high ratio gear set shown in 
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Fig. 9. Variation of FE Ratio with Face Width. 


Example 3, a considerable percentage variation can exist in Lmin 
depending upon the pitch selection. These comparisons demonstrate 
the importance of careful consideration for critical designs. 


3. Stress Pattern 


Bending Stress 


A most significant difference exists between the stress pattern 
developed in spur and helical gear teeth because of the oblique line - 
of-contact. In the spur tooth with uniform load distribution the 
structure can be normally considered as a beam and the classical 
Lewis analysis applies. 

With the oblique line-of-contact any approach using a beam 
analysis can result in serious error. The correct theory to use 
is that of a cantilever plate. According to the principles of a can- 
tilever plate for any position of load the value of the maximum 
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Table 
Examples of L.,;, Variations 
Numbers | 
Diametral | of Teeth ; 
| Example Data Pitch (Pinion) | Lin” Remarks 
m, - 2.0 1 37 | 14.8 | Changes in L,,;, 
1 contien cam | are relatively small 
ae and the influence of 
face - 10” 2 74 | 16.1 pitch or number of 
| teeth in the pinion 
3 | ill 15.8 is minor 
+ = 4 + 
m, - 2.0 1 10 | 12.5 Low numbers of 
| | 
~ pinion teeth with 
o | 
. comtare - 15.8 coarse pitches lead 
| face - 10” | 2 20 | 13.9 to a considerable 
| percentage varia- 
3 30 14.6 tion in total contact 
4 40 15.1 length 
+ + , 
m, - 10 1 | 10 | 12.8 Even with high 
_ 552° | ratios L__;,, varies 
3 | centers - 554 | | | qansttesell do. 
face - 10” 2 | 20 14.4 pending upon pitch 
| 3 30 15.0 
| | 4 | 40 17.1 











moment at the fixed edge is independent of the absolute width of the 
plate or the absolute normal distance between the load and the fixed 
edge, being dependent only on the relative value of these dimen- 
sions. This is explained in the paper by Seireg and Wellauer (1).' 
The effect of the inclination angle on the moment distribution at the 
root is illustrated in Fig. 10. 

It is evident that the “worst” position of loading for maximum 
stress can vary depending upon numerous conditions. Fig. 11 shows 
the moment distribution for various positions of the load line. The 
effect of buttressing and of change in load due to tooth end deflection 
or contact relief can also modify both the magnitude and the position 
of the maximum stress. 


Contact Stress 





The oblique line-of-contact also manifests its nature in regard 
to the contact stresses which determine the surface durability or 


‘Numbers in parentheses designate References at end of paper. 











12 E. J. WELLAUER HE 











1.0 
MOMENT = Cah | Tt 
@ «= LOAD INTENSITY co 
ON CONTACT LINE 
8 h = TOOTH HEIGHT po 
cu 


. Pt 
/ 
/ 





- FACTOR 








c 

om 

—— 
Sc 























ta i 


h 


i 


20° 











Fig. 10. Effect of Inclination Angle on the Moment 
Distribution at the Root (Load Line Away from the Corner). 


resistance against pitting. Generally, spur gears can be considered 
as cylinders in contact whose relative diameters vary with the po- 
sition of the contact along the profiles. The helical gear on the 
other hand must be considered as two truncated cones in contact. 
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Fig. 11. Helical Gear Contact Conditions and Bending Moment. 
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This is a rather complex stress system which, as yet, has not been 
completely solved. The worst position for spur gears will be at that 
position of single tooth contact having the smallest relative radii of 
curvatures. On the other hand, for helical gears it is obvious that 
varying relative radii of contact exist simultaneously and experience 
alone has determined which value is important; this particularly be- 
cause of the significant influence of elasticity on the distribution of 
the load. 

The face of spur and helical gears are often crowned for better 
operation. The extension of the localized contact is proportionate to 
the cube root of the load for both types of gears. 


4. Load Distribution 


The load distribution between teeth and along the contact lines 
is very important. The oblique line-of-contact by virtue of its 
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Fig. 12. Helical Gear Load Distribution Factor -K 
















































































Fig. 13. Effect of Different Types of Loading 
the Bending Moment Distribution at the Root 
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contact action and uniformity of total length importantly affects the 
load distribution. 

An effect is present on tooth deflection, the distribution of load 
along the individual lines-of-contact, the thermal expansion and the 
difference between cylindrical and conical contact. The analysis of 
the effect of tooth deflection is presented by Wellauer (2) and is pri- 
marily concerned with the fact that for spur gears the maximum 
rigidity affecting tooth contact is normally secured with two teeth in 
mesh. For helical gears, contact conditions are more uniform and 
the distribution of loading between two teeth is dependent upon the 
ratio of the length of action in the transverse plane to the base 
pitch. This is illustrated in Fig. 12. 

Deliberate profile and lead modifications to affect the load dis- 
tribution across an individual line-of-contact can result in varia- 
tions of bending moment illustrated by Fig. 13. 
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5. Elastic Characteristics 


The special elastic characteristics pertinent to the oblique line - 
of-contact of helical gears has practical significance both in regard 
to the performance of the blank and of the tooth. These phases in- 
volve vibration, dynamic loading and the width and depth of the 
flattened area. 

Analyses of load sharing, deflection of teeth, and blank deflec- 
tions are available in the literature. 

Insofar as vibration is concerned, helical gears have been used 
almost universally for high pitch line velocities. There are several 
reasons for this preference. First, the smooth action previously 
discussed minimizes tooth excited forces and noise. Secondly, since 
the contact conditions are more uniform, there is less excitation 
due to the change of elastic mesh stiffness such as occurs with the 
one to two tooth load cycle in spur gears. 

The oblique line-of-contact has an influence upon the dynamic 
loading developed by the gear teeth themselves. Fig. 14 based upon 
Harris (3) shows the derating expressed as a dynamic factor caused 
by the effects of base pitch errors and the load with reference to the 
pitch line velocity. Aiso plotted are commonly used velocity factors 
applied for rating commercial gears. Curves 1 and 2 are used for 
helical gears and precision spur gears whereas 3 and 4 have been 
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Variation of Dynamic Factor for Spur and Helical Gears. 
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used for less accurate spur gears. These curves follow the trend 
which might be predicted by the analyses. 

It also is obvious that because of the helical shape, the metal in 
a helical gear blank is distributed more symmetrically around the 
cylindrical surface and opportunities for exciting lateral vibrations 
are at a minimum. 

Fig. 15 is an example of the elastic contact deflections which 
can be obtained across the line-of-contact of helical gears. It is to 
be noted that the flattened land is of a significant width even for ex- 
tremely minute deflections normal to the profile. 
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Fig. 15. Profile Deformation Along Contact Line. 


Conclusion 


Beside the factors discussed in this article others are equally 
important and sometimes deserve detailed analysis. The nature of 
the helical gear oblique contact line is a controllable parameter 
which can be usefully employed for producing optimum designs. 
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Basic Geometry and Tooth Contact of 
Hypoid Gears 


By Meriwether L. Baxter, Jr. 
Chief Research Engineer 
The Gleason Works 
Rochester, New York 


Basic Geometry of Hypoid Gears 


There is no mystery about hypoid gears. There can hardly be 
any real mystery about a machine element that is produced in the 
millions each year by many unrelated companies; hypoid rear 
axle gears perform with quiet competence in vehicles all over the 
world. 

The series of articles by Ernest Wildhaber in American Ma- 
chinist in 1946 contains a complete description of the geometry 
and kinematics of hypoid gears; since the publication of the Wild- 
haber articles there has been acertain amount of correspondence 
with engineers in Europe and Japan indicating that they have made 
an intelligent study of at least parts of the series. There have 
been relatively few similar inquiries from American engineers, 
which perhaps indicates that the presentation was somewhat too 
general and too complete to appeal to the engineer with a production 
job to do. 

This paper presents a simplified but logically complete de- 
scription of hypoid geometry and kinematics which should clarify 
the meaning of the various terms encountered in hypoid design and 
calculation, and the relationships involved. 


Synthesis of Hypoid Geometry 





Our geometry contains at the beginning two axes in space; in 
general they do not intersect, and are neither parallel nor or- 
thogonal. Now we can define a new line which is the common per- 
pendicular to the two axes, called the line of centers. The distance 
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Figure l. 


between the two axes measured along the line of centers is the 
offset E. Figure 1 shows the axes X, and Xp and the line of 
centers O,Q,. The angle between the axes is naturally measured 
in a plane parallel to both, and its symbol is £. For definiteness 
we shall call X, the gear axis and Xp the pinion axis. 

We now propose to transmit motion from one axis to the other 
by direct tangent contact between a pair of surfaces, that is, by 
gear teeth; and it is clearly necessary to make a decision as to 
where in our space the contact is to take place. This decision is a 
critical one, and determines the kind of gears we shall end up with. 
Special cases yield special results; for example if the tooth mesh is 
to take place on the line of centers we shall obtain worm gears. 
Since hypoid gears are the most general kind of gear, we must 
choose a general point, such as M in Figure 1. 

The mean point M can be specified in a variety of ways; we 
choose to locate it by the three quantities, R, ¢, and Zp shown in the 
figure. It will be seen that Zp specifies the plane of rotation of the 
gear in which M lies; that R is the mean radius of the gear; and 
that ¢« gives the orientation of M in its plane of rotation. 

Through M we can construct line CD intersecting both axes; this 
line is called the pitch vertical and is an important element in our 
geometry. It will be seen that the inclination of the pitch vertical 
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Figure 2. 


when viewed along the gear axis Xc is the given angle ¢€; similarly 
the angle 7 is the inclination of the pitch vertical when viewed along 
the pinion axis. 

Rp is the distance of mean point M from the pinion axis; that is, 
the mean pinion radius. In the material that follows we shall as- 
sume the axes to be at right angles in order to simplify formulas 
and diagrams; the principles however are general. 

From the geometry of Figure 1 the formulas for the newly de- 
fined quantities 7 and Rp are as follows: 


tan) = = (B - Reine) (1) 
Zp 
Rp = Zp sec 7 (2) 


Let us now consider a plane containing mean point M and per- 
pendicular to the pitch vertical. Figure 2 shows such a plane, inter- 
secting the gear axis at H, , the gear apex, and the pinion axis at H,, 
the pinion apex. For reasons to be seen later we shall call this the 
bitch plane, the distance MH, the gear cone distance A, and the dis- 
tance MH, the pinion cone distance Ap. Likewise the angle between 
the gear axis and the pitch plane is the gear pitch angle IT, while 
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the pinion axis is inclined at the pinion pitch angle y to the pitch 
plane. Other quantities to be defined are the angle « = { H, MH, in 
the pitch plane and the linear dimensions Z and G locating the gear 
and pinion apices respectively. 

From the geometry it is possible to write many relationships 


between our defined quantities; those listed below are sufficient and 
relatively simple 


tany = sin 7 ctn « (3) 
cos [= siny/cose (4) 
Sine’ = sin e/cos ? (5) 
A= R/sinI (6) 

Ap = R,/sin? (7) 
Z= Acosr -Zp (8) 

= -Rcos€ + Ap cos? (9) 


Since the pitch vertical intersects both axes, and the pitch plane 
is perpendicular to the pitch vertical, it is possible to describe a 
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pair of cones, one about each axis, which are tangent to each other 
at M. The apex of the gear cone is at H, and that of the pinion cone 
is at H,; the lengths of the conical elements to the mean point M are 
the cone distances A and Ap previously defined. Figure 3 shows 
these pitch cones. 

It should be emphasized that these pitch cones are not “pitch 
surfaces” in the sense that this term is used for gearing with paral- 
lel or intersecting axes. For example, we cannot say that the pitch 
cones roll together without sliding, since up to this point we have 
not introduced the velocity ratio. The pitch cones should be con- 
sidered merely as geometric elements to be used as the basis of 
design of the gear blanks. 

To summarize to this point, we have started with a pair of axes 
related by the quantities E, the offset, and £, the shaft angle. We 
then assumed the location of a mean point M of the tooth mesh, de- 
fined by the quantities R,¢«, Zp. Through the mean point we have 
defined a pitch plane, and two pitch cones tangent to each other and 
to the pitch plane. We have determined formulas for the pertinent 
quantities, including the location of the apices, the cone distance, 
the pitch angles, and the angle e’ in the pitch plane between the two 
cone elements. 


Velocity Ratio Determines Spiral Angles 





Figure 4 is a view of the pitch plane; we are looking along the 
pitch vertical from the pinion axis toward the gear axis; H, and H, 
are the gear and pinion apices respectively. The gear pitch cone 
lies below the view, and the pinion pitch cone above. The vector ¢g 
is a unit vector in the direction of the gear axis; associated with it 
is the gear angular velocity “Gc. Similarly vector p has the direc- 
tion of the pinion axis and is associated with the angular velocity w,,. 
Vectors A and A, run from the apices to the mean point M; their 
magnitudes are the cone distances A and A ps 

The vector velocity of M considered as a point on the gear is 


Vo =% gxA (10) 
and its magnitude is A sin '-,.. Similarly for the pinion 
Vp = #, pxA, (11) 


of magnitude A, siny-,. Due to the assumption of direction 
chosen for p and g, the relative velocity vector may be written 


pxA, (12) 
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The relative velocity vector V lies in the pitch plane, and serves 
to determine the tooth direction. If tooth contact is to take place 
at M with the specified velocity ratio the fundamental law of gear 
tooth contact must hold: the relative normal velocity must be zero. 
This law is expressed very simply in vector notation 


V-n=0 (13) 


where V is the relative velocity vector and n is a vector along the 
common contact normal. 

As seen in Figure 4, we define the gear spiral angle y as the 
angle between the tooth tangent and the gear pitch cone element, and 
the pinion spiral angle Vp as the angle between the tooth tangent and 
the pinion pitch cone element. 

The sliding velocity for contact at M is the magnitude of vector 
V and can be calculated by 


v.= V¥.V (14) 
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The spiral angles may be obtained from 
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cosw = . . , (15) 
s 

cosy, = te < (16) 
Ss 


The more familiar algebraic formulas are 





cosy _A siny N_R,_N 

cosy, AsinE “n Rn (17) 
Vp=vtre’ (18) 
V; = Rwo cos y (tan Wp - tany) (19) 


The tooth normal may have any direction perpendicular to V or 
the tooth tangent; thus the pressure angle is subject to free choice. 
The spiral angles however are fixed by the chosen mean point and 
velocity ratio. If the spiral angles are not satisfactory we must 
change one of the assumptions for the location of the point M; usu- 
ally Zp is changed for this purpose. 


Surfaces of Action 





Within reason we can choose any tooth shape for one member, 
say the gear, and the conjugate surface for the mate and the surface 
of action can be determined. If the chosen mean point M is to lie 
in the surface of action the direction of the tooth tangent must of 
course be as calculated above. 

The suitability of a particular tooth shape can be studied by 
reference to its surface of action. The relationship between the path 
of action and such conditions as undercut and contact ratio is well 
understood, and it is natural that similar considerations should 
govern the choice of tooth surfaces and surfaces of action for hypoid 
gears. In particular we may wish to specify (a) that the two sur- 
faces of action be equally inclined to the pitch plane, and (b) that the 
two surfaces of action intersect in a line tangent to the pitch plane. 

The first of these aims is intended to produce equal contact ratio 
and equal freedom from undercut for the two sides of the teeth. For 
good reason we sometimes wish to depart from this balanced condi- 
tion, but to do this intelligently we must first know how to get the 
balance. Failure to meet the second requirement at least approxi- 
mately will result in teeth which are pointed at one end and undercut 
at the other. 
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Figure 5. 


Due to the shaft offset of hypoid gears an unsymmetrical condi- 
tion exists such that the surface of action makes an angle with the 
pitch plane which is not equal to the pressure angle of the tooth sur - 
face. In fact there is a particular hypoid pressure angle, called the 
“limit pressure angle” which produces a line in the surface of action 
tangent to the pitch plane, and which therefore may be said to be 
equivalent to zero pressure angle for other types of gears. 

Figure 5 is similar to Figure 4, and shows a unit normal vector 
to the gear tooth, f,, inclined to the pitch plane at the limit pres- 
sure angle 6,. As the gear rotates through a small angle @, point M 
of the gear is carried to M’, where 


MM’ = SA = 6gxA (20) 
SA = OA (21) 
The change in n, due to rotation @ is 


An =6@¢xn, (22) 


If M’ is to remain a point of contact, then the formula V - n = 0 must 
continue to be fulfilled, or 


V- An+AV-n,=0 (23) 


Differentiating equation (12) we obtain 


AV = 4 gxdA+w px SA, (24) 
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Substituting from equations (12), (20), (21), (22), (24) in equation (23) 
gives 


wo gx A+w, px A,)-(@ gx No) + WG E+ Wp p) x (6 Ex A)-No =0 


which can be simplified to 


no - [px (g@xA) - gx (PxA,] = 0 (25) 


This is sufficient to define the limit pressure angle ¢,. 
The equivalent algebraic formula for ¢, is 


_ Ap sin Wp - A sin 


~ RG * Aptany + Atanr (26) 





This value of 6, is regularly used in hypoid design to determine 
the amount of unbalance for the gear tooth pressure angles. 

Now let us assume a gear tocth of the limit pressure angle 6, to 
be a surface of revolution about an axis perpendicular to the pitch 
plane, as seen in Figure 6. Unit vector c is directed along the axis 
of revolution, and a vector (rp,.n,) connects the axis to mean point 








an 6é 1 o "se 
af Ge « A, 
Pr we 6 * we B) * tae 8 
o 6 - x x Ne 
4 a, an 
2 
ve 
ne + 
fre io «x V 
On We o ow 
o COS o,) — —— ne —T 
zo on $ |——— ony | 
Ap COS Wp A cos w \Ap tan y A ton T/ 











28 MERIWETHER L. BAXTER, JR. 


M. If we rotate M about c to a new point M”, the latter point re- 
mains on the tooth surface and in the pitch plane. If the angle of 
rotation is 6, 





BA = 6cxr,,n, (27) 
An = @cxn, (28) 
as before, AV = (wg + wp) x AA (29) 
and AV -n, + V- An, = 0 (30) 


By substitution, as above 


(WG & + &, P)xX (6 tno CXMo)- No + V-Gcxn, = 0 (31) 


If we let (wc g + w, p) = Ww =the relative angular velocity, then this 
relation reduces to 


V2 


= eS 
Ino - Ss: (a x V) (32) 


Thus there is a particular normal radius of curvature of the 
gear tooth which when used with the limit pressure angle produces 
a surface of action tangent to the pitch plane. Since @, is in effect 
an average pressure angle between the two sides, the radius of 
curvature thus calculated is a good choice for real teeth of ordinary 
pressure angle. 

The usual algebraic formula for the limit radius of curvature is 


tan Vp - tany 
(no COS b,) = nent (33) 
1 1 "' tany, tany 


et a en 
Po A, tan? ” Atanr’ 











A, cosy, Acosy 


In order to arrive at a limit radius of curvature corresponding 
to a standard size of cutter used for manufacturing hypoid gears, we 
can vary the assumed value of angle ¢«, which was one of the param- 
eters originally used to specify the location of the mean point M, 
just as we varied the assumed value of Zp to obtain suitable spiral 
angles. While this trial-and-error procedure sounds cumbersome 
it is not difficult in practice, and a satisfactory design can be 
achieved in a short time with the use of a desk calculator. 
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It is of course possible to calculate actual points in the surface 
of action for a particular tooth surface. To demonstrate the effec- 
tiveness of the concepts of limit pressure angle and limit radius of 
curvature, we have calculated these values for a typical automotive 
hypoid ratio; we have chosen pressure angles unbalanced by exactly 
the limit pressure angle, and radii of curvature exactly equal to the 
limit value. 

Figure 7 shows the trace of the surfaces of action in both the 
pitch plane and the normal plane. It will be seen that the surfaces 
of action are only moderately curved, and that the pitch-plane traces 
lie close together, indicating that the intersection of the surfaces of 
action lies very close to the pitch plane. In the normal plane we see 
that the lines of action are well balanced for direction. 

These results justify the use of the limit pressure angle and the 
limit radius of curvature as basic design parameters for ordinary 
cases. In extreme cases of offset, ratio, or spiral angle we do of 
course make a detailed study of the actual surfaces of action. 
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Mismatched Gears 





A primary requirement of real gears is that they perform well 
under actual operating conditions. Owing to manufacturing toler- 
ances and deflections under load, gear and pinion are seldom if ever 
operating in their theoretical relative positions; nevertheless we 
expect them to transmit the required power smoothly and quietly. 

Theoretically conjugate bevel and hypoid gears can be manu- 
factured, but they are not practical because of their sensitivity to 
the misalignments inherent in any real gear housing. In order to 
obtain the required adjustability we introduce into the teeth of one 
or both members carefully chosen modifications usually referred to 
as “mismatch.” By properly choosing the kind and amount of mis- 
match best suited to a particular application, we obtain results far 
superior to those which would result from theoretically conjugate 
teeth. 


Tooth Bearing 





It is a general practice to judge the type and amount of mismatch 
in a pair of gears by running them together in a testing machine 
under light load with the tooth surfaces painted with a marking com- 
pound such as red lead. The resultant marking on the tooth is called 
the “tooth bearing.” The shaded area in Figure 8 represents a tooth 
bearing. 

Correctly made teeth will show a tooth bearing centrally located 
on the tooth surface, when the gears are run in their standard rela- 
tive positions. The amount of lengthwise mismatch is readily ob- 
served by the length of the tooth bearing; under light load the bear- 
ing will frequently cover about one half of the face width. Profile 
mismatch is not observed as readily owing to the fact that the tooth 
height is relatively small; nevertheless the existence of profile 
mismatch is evidenced by the lack of any heavy concentration of 
bearing at the tips of the teeth. 

By displacing the pinion specified amounts on the testing ma- 
chine we observe the resultant movement of the tooth bearing, and 
are then able to judge the adequacy of the set to operate in its par- 
ticular application. “Development” of a gear pair is the process of 
generating pinions with a variety of combinations of mismatch to 
determine the one best suited to the job. 

The testing machine adjustments which have been usually used 
in the past to displace the bearing are a change in the running offset 
and a change in the pinion axial position. The former is frequently 
called a “vertical” change and represented symbolically by V. The 
latter is called “horizontal” change and represented by H. The 
combined V-and-H change required to cause the tooth bearing to 
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Figure 8. Figure 9. 


move from one end of the teeth to the other while remaining at an 
average height on the profile is called the “V-and-H check” and is 
frequently used to characterize the over-all nature of a particular 
tooth contact. This is frequently written in the form of a fraction; 
for example 


means that a combined change of 0.021 in. vertical adjustment and 
0.015 in. of horizontal adjustment on the testing machine will move 
the tooth bearing from toe to heel, while maintaining a good profile 
position. However, these figures do mot indicate the amount of in- 
dividual vertical or horizontal adjustment which can be tolerated. 

The tooth bearing also serves to detect the existence of a con- 
tact condition known as “bias” in which the tooth bearing extends in 
a diagonal direction across the tooth. Some degree of bias is some- 
times introduced to obtain a quieter job, or to produce a different 
V/H ratio. Like lengthwise and profile mismatch, bias can be read- 
ily controlled by alterations in the setting used to generate the 
pinion. 

Newer designs of testing machines use a change in shaft angle 
instead of the vertical to accomplish a shift in bearing from end to 
end. 

By listening to the gears in the testing machine, observing the 
position and shape of the tooth bearing, and observing the direction 
and amount of tooth-bearing movement under displacement, it is 
possible to develop a gear pair tailored to the particular require- 
ments. The method in general is not however entirely satisfactory 
for the following reasons: 


1. It is relatively costly and time-consuming. 
2. Personal judgment is an important factor. 
3. It provides little basis for further investigation. 
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For these reasons we have developed an analytical method for 
determining quantitatively the conditions of contact and motion 
transmission in mismatched gears; this has been called TCA (tooth 
contact analysis). 


Brief Description of TCA 





TCA permits us to calculate the tooth-bearing characteristics 
for a bevel or hypoid gear pair in any specified running position; in 
addition, it provides graphs of the transmitted motion. Thus to 
some degree it replaced conventional development using the testing 
machine. 

On the other hand it is entirely theoretical, and does not at 
present consider such important practical matters as bending and 
compression of teeth under load, or distortion due to heat-treatment. 
For this reason TCA should be considered a supplement to and not a 
replacement of traditional methods. An important advantage how- 
ever is that it provides a sound theoretical basis for further re- 
search in the subject of gear tooth action now being carried on at 
the author’s company. 

The data required for TCA are simply the machine settings used 
to manufacture gear and pinion, and the specified running position of 
the pair (not necessarily the nominal running position). 

The results of TCA are as follows: 


1. Data for plotting the tooth bearing. 
2. Data for plotting motion graphs. 
3. Calculated V-and-H check. 


Calculated results of TCA correlate well with observed results 
under light loads in the testing machine; moreover, rather small 
changes which would be difficult for a machine operator to set or 
observe produce definite and repeatable effects. 


Qualitative Description of Tooth Contact 


In order to establish a theoreticai basis for the study of tooth 
contact it is first necessary to arrive at a description of what actu- 
ally takes place as gear teeth mesh together with combined rolling 
and sliding action. To this end we first discuss conjugate gears, 
and then show the effects of mismatch. 


Fully Matched Gears 


Let us consider real gears rotating about fixed axes with a con- 
stant velocity ratio, and equally spaced identical teeth. 
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We may take any reasonable surface as the tooth surface for one 
member and determine the conjugate surface for the other member, 
such that (a) the velocity ratio is maintained, (b) there is line con- 
tact within the physical boundaries of the tooth surfaces. 

As rotation proceeds, the line of contact progresses across the 
tooth surfaces of both members, usually changing direction and 
shape as it does so. The angle of rotation during which a particular 
line of contact exists is determined by the tooth boundaries, and is 
greater than one pitch in properly designed gears. The ratio of this 
angle to the angular pitch is known as the “contact ratio.” 

Figure 9 shows a tooth surface and the successive positions of 
the line of contact. The solid lines represent two positions sepa- 
rated by a rotation of one pitch. The ratio of the number of lines 
encompassed by the tooth surface to the number for one pitch is the 
total contact ratio. 

Figure 10 is a motion graph for conjugate teeth. The ordinate 
represents error in position of the driven member, zero in this 
case. The abscissa represents rotation of the driver. It will be 
seen that the action of successive teeth overlap and that there are 
two teeth in contact during the shaded intervals. This can occur 
because the velocity ratio is constant. 


Single Mismatch 


Now consider the tooth surface in Figure 11. This surface was 
obtained by starting with the tooth in Figure 9 (which was conjugate 
to its mate) and relieving it on both sides of line AB in such a way 
that the surface remains continuous and tangent to the original sur- 
face along line AB. 

As this tooth meshes with the mating member, there will bea 
position where the line of contact with the conjugate member was 
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Figure 10. Figure ll. 
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at (1). Now, however, only point Al can be in contact, since all other 
points on the line of contact have been relieved. Since point Al is a 
point of the original surface, the position of the driven member 
relative to the driver is correct. Therefore: 


1. There will be point (not line) contact. 
2. Uniform velocity will be maintained. 
3. The path of the contact point is line AB. 


Figure 10 is also a correct motion graph for single mismatch. 


Combined Mismatch 





Most spiral bevel and hypoid gears are relieved both profilewise 
and lengthwise, so that only a single point in the center of the tocth 
is unchanged from the theoretically conjugate surface. Combined 
mismatch permits the gear pair to tolerate displacements under 
load and assembly errors. 

With combined mismatch, as with single mismatch, there is the- 
oretically point contact at any instant, and a path of contact across 
the teeth. However, the velocity ratio will not remain constant, and 
there is a theoretical instantaneous velocity change at each tooth 
changeover which is of interest to us analytically although its actual 
effect will be largely damped-out under load. 

It is the purpose of TCA to reveal the contact conditions for 
combined-mismatch teeth. 


Motion Curves for Combined Mismatch 





The upper part in Figure 12 is the same kind of graph as Fig- 
ure 10. We have plotted the error in displacement of the driven 
member against rotation of the driver for the entire duration of 
contact for one pair of teeth, as though no other teeth existed. 

As would be expected, the driven member lags behind its cor- 
rect position except at one instant. 

The lower part in Figure 12 is the corresponding velocity curve. 
The velocity of the driven member is greater than standard at the 
start and less at the end. 

The curves in Figure 12 are typicai of teeth with profile mis- 
match. 

Figure 13 shows the motion curves in Figure 12 with successive 
tooth meshes superimposed. 

The intersection of the successive displacement curves is the 
“transfer point” at which contact shifts from one tooth to the next. 
The transfer points are one pitch apart. 
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The portions of the displacement curves beyond the transfer 
point are therefore not used, but their length is an indication of the 
amount of extra contact area available on the teeth to provide ad- 
justability and accommodate misalignment. 

The displacement error at the transfer point is the theoretical 
amount the driven member oscillates with each tooth mesh. This is 
usually a few tenths of a thousandths of an inch. 

The velocity curves in Figure 13 do not intersect; at the transfer 
points there is a sudden jump in the velocity of the driven member, 
which is the velocity with which the new driving tooth strikes its 
mate. This value is therefore a measure of the impact occurring at 
each tooth mesh. 

As profile mismatch is increased, the impact and oscillations 
increase in magnitude, resulting in noise and stress. It is therefore 
essential that mismatch be kept to the smallest value compatible 
with the adjustability requirements. 

It should be emphasised again that these motion errors are only 
theoretical concepts useful for analysis and research, and that the 
real motion transmitted will be essentially uniform due to the aver- 
aging effect produced by the application of a load, and its resulting 
deflections. 








Contact Mark and Tooth Bearing for Combined Mismatch 





Theoretical point contact is, of course, only an imaginary con- 
cept. Under even light load, tooth surfaces compress and deflect to 
cause surface contact. Furthermore, the grain size of any marking 
compound used is sufficient to show area contact around the theo- 
retical point of contact. 

It has been found that gears on a testing machine under light 
load will show marking wherever the surfaces are closer than about 
0.00025 in. (one quarter of a thousandth). Thus we can calculate the 
shape and size of mark which will show around a particular point of 
contact. 

These contact marks are roughly elliptical, with the long axis in 
about the direction of the lines of contact of the conjugate gears. 
Their length is determined by the amount of mismatch. Such marks 
can be observed easily on the testing machine with a suitable mark- 
ing compound by holding one member stationary and rocking the 
mate against it. 

As gears roll together in the testing machine, the contact marks 
move across the tooth surfaces to produce a total marking known as 
the tooth bearing. 

Figure 14 shows how we may indicate the calculated bearing. 
Points Tl and T2 are points of contact corresponding to transfer, 
and the path of contact is indicated as a slightly curved line between 
Tl and T2. 

The elliptical contact marks have been indicated for points T1 
and T2 and for another contact point in between. 

The theoretical tooth bearing, therefore, is bounded as follows: 


= portion of contact ellipse for T2. 

= deepest contact of teeth of mate. 
envelope of ends of contact ellipses. 
= portion of contact ellipse for T1. 

= tip of tooth. 

= envelope of ends of contact ellipses. 
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In the case shown in Figure 14, the entire path of action from 
transfer point to transfer point lies on the tooth surface. There- 
fore, we do not expect the marking at E and B to be particularly 
heavy. 

In Figure 15, we have a case where one transfer point, T1, lies 
outside of the tooth boundaries. As contact moves from T2 to X, 
conditions are similar to those in Figure 14. Beyond X, however, 
motion must be transmitted by the edge of the tooth E to the left 
of X, resulting in a heavy edge contact line and a distortion of the 
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Figure 14, Figure 15. 


bearing in this area. This is the sort of effect produced by an axial 
displacement of the pinion too large for the amount of profile mis- 
match provided. 

Figure 16 shows the form in which TCA results are usually pre- 
sented. The analysis has been performed for three separate running 
positions corresponding to toe, central, and heel contact as indicated 
by the three V/H values shown. Displacement curves and theoreti- 
cal tooth bearings are shown for each running position. 

This particular diagram shows practically perfect uniformity of 
motion transmission, since the scale used to plot the curves is very 
large; the value of 0.00002 radians indicated on the ordinate scale 
represents less than 0.0001 in. in peripheral motion of the gear. 
There is overlap in all cases, but little adjustability since a small 
amount of axial adjustment will change the slope of the curves and 
result in tip contact. 

The tooth bearings are practically perfect; the paths of contact 
are straight and vertical, and the length and inclination of the con- 
tact marks are constant. 


Mathematical Basis of TCA 





A fundamental requirement of TCA is that the given tooth sur- 
faces be generated in a completely arbitrary manner; that is, we 
must be prepared to analyze a pair of gears made with machine set- 
tings arrived at by some unknown procedure, or by trial-and-error. 
Similarly we must be able to analyze operation in relative running 
positions different from those specified on the drawings. 

These requirements imply that we cannot make the assumption 
of uniform velocity transmission which so simplifies ordinary gear 
mesh analysis, and forces us to adopt a completely general proce- 
dure. As this procedure is described it will become evident that the 
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Figure 16, 


calculation would be hopelessly long without the use of a high-speed 
electronic computer. 

In general a point on a tooth surface can be specified by two pa- 
rameters. In the case of generated tooth surfaces these parameters 
are taken to be the cradle angle q and an angle specifying rotational 
position around the cutter axis, ¢; for FORMATE ® or HELIXFORM® 
teeth convenient parameters are the distance s from the tip of the 
cutter blade to the specified point, and an angle a about the cutter 
axis. 

If we assume a pair of values (q, 6) we can determine from the 
conditions of generation all we need to know abouta specific point on 
the pinion tooth, including the direction of the tooth surface normal 
at that point. Similarly the choice of a pair of values (s, a) locates 
a point on the gear tooth with its normal. Since these points have 
been chosen arbitrarily, we have no assurance that they would ever 
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come into contact with the gears running in the specified relative 
position. 

Nevertheless we can imagine the two gears to be physically 
moved into such a position that the two chosen points and their 
normals are coincident. Figure 17 indicates such an imaginary 
juxtaposition. Furthermore, we can imagine one of the gears to be 
rotated about the common surface normal, maintaining contact at 
the chosen points, until the angle between the shafts becomes the 
specified shaft angle >. 

When this has been done the gear and pinion blanks are in afixed 
relative position, and we can calculate the shaft offset E, and the 
axial positions xp and xq. These calculated values will not of 
course be the specified running positions, so we now proceed to 
change the assumed values of the three parameters 6, a, s until the 
values of E, xp and xq are correct, and thus obtain a point of con- 
tact. At this point of contact we can work out such details as the 
angular position of the gear and pinion, the instantaneous velocity 
ratio, and the curvature characteristics of the contacting tooth sur- 
face, that are required as part of the results of TCA. 

Now by repeating the procedure with a different value of q we 
can obtain the sequence of contact points which makes up the path of 
contact. 





Figure 17. 
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Figure 18. 


Figure 18 will serve to clarify the sequence of operations just 
described. It will be seen that we require some iterative method to 
vary three parameters until three calculated errors, or residuals, 
become zero to some specified limit. The next section describes 
such a procedure. 


Multiple Iteration Method 


It is not particularly difficult to handle multiple interpolation 
when the relating functions are simple and short enough to make 
differentiation easy. In many engineering problems however this is 
not the case and we need a method which does not require any analy - 
sis of the functions involved, but merely the ability to calculate them 
rapidly. 

Let us say that if we are given values of three parameters a, b, 
c, we can calculate three quantities, x, y, z which depend on a, b, c. 
It is not however possible to calculate a, b, c from given values of 
x, y, Zz. Nevertheless we must find somehow what values of a, b, c 
to choose to make x = y = z = 0, and we know as a practical fact that 
such a solution exists. 

If we make a small change in a and determine its effect on x, y 
and z, we have obtained approximations of three partial derivatives 
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Similarly by stepping b and c small amounts we can complete the 
matrix of 9 partial derivatives or rates, and with these we can cal- 
culate new values for a, b, c merely by solving three simultaneous 
linear equations. Repeated application of this procedure will pro- 
duce convergence, if 


: a) there is a unique answer 

b) the starting values of the parameters are reasonable 

c) high accuracy is maintained in computation, particularly in 
the solution of the linear equations. 


Figure 19 is a computation flow-chart for one way of handling 
such an iterative procedure; not all necessary details are included, 
but only enough to indicate the plan. 

We begin by initializing the starting values of a, b, c and speci- 
fying the small steps Aa, Ab, Ac to be used to get rates in the first 
pass. We also set switch 1 for the first pass. 

The computational block calculates x, y and z, and the test block 
determines whether the errors are sufficiently small. If so, we exit 
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from the iteration program. Usually the test will fail and control 
will pass to switch 1. 

Branch (1) of switch 1 resets itself to branch (2), proceeds to 
set switch 2, step a, store last answers and return to computation. 
On subsequent passes a new set of rates is figured following each 
change in a, b, or c. 

When three sets of rates have been determined, new values of 
Aa, Ab, Ac are calculated which are first approximations to a solu- 
tion. If we had reason to believe the functions are linear it would be 
desirable to substitute all of these changes at once; however it is 
generally more practical to apply the changes one at a time as was 
done with the original steps, obtaining new rates at each stage in 
preparation for subsequent passes. 

This procedure, for one, two, or three variables, has been of 
tremendous value in our engineering computation; however it must 
of course be applied intelligently. One particularly troublesome 
condition is unsuspected discontinuity in the functional relationships; 
if such discontinuities exist they may be located in advance, and the 
program modified to prevent stepping across this region. 

Multiple solutions have not been troublesome since their exist- 
ence can usually be predicted, and sufficiently good starting vaiues 
will force the solution we are after. 

When the iteration does not converge, and the program has been 
carefully written to avoid round-off errors, we usually find that 
there is a physical meaning corresponding to the non-convergence. 
For example, in the TCA program, failure to converge indicates 
either 


a) line contact instead of point contact, or 
b) contact at or near the cusp of a generated surface. 


Conclusion 


This paper has developed a geometrical frame of reference for 
the study of hypoid gears; from this basis we have proceeded to a 
discussion of tooth contact in conjugate gear pairs with line contact, 
and then to the more important practical case of tooth contact in 
gear pairs so modified as to permit motion transmission under con- 
ditions of misalignment. As a necessary tool in the discussion, an 
iterative procedure has been described which has wide applicability 
in engineering work. 

We hope that the material discussed here will have theoretical 
and practical value to both mathematicians and engineers. 








Stress Concentration Caused by a Small 
Circular Hole in a Bent Plate 


By Lajos Imre Nagy 
General Motors Corporation 


Consider a flat plate (see Fig. 1) bent at the ends by the mo- 
ments (M). A concentric hole at the center of the plate will effect 
the stress distribution in the neighborhood of the hole by causing 
stress concentration. According to Saint-Venant’s principle the 
effect of a small hole may be neglected at a distance of a few di- 
ameters from its edge. The points which are at such a distance 
from its edge can be regarded as being at infinity. 

In order to solve the two-dimensional problem for finding the 
stresses at any point in the plate the equation of equilibrium, the 
compatibility equation,* and the boundary conditions must be satis- 
fied. 

By introducing the stress or Airy’s Function y, the equations or 
equilibrium are automatically satisfied by taking the following ex- 
pressions for stress components: 














Two-dimensional cartesian Two-dimensional polar 
coordinate system coordinates 
: 2 “ae ret - oe 
= ay as ror r2 aeg2 ~’ 
. - 0? ee _- _0*y 
iggy Be Tate 9 % b, >(1) 
~ = os de 2 (124) ¥ 
xy ~ @xay ‘r@ ~ @r\r 06 —~ 





By substituting the stress components into the compatibility 
equation, one may find that the stress function must satisfy the bi- 
harmonic equation: 

" *§, Timoshenko and J. N. Goodier, “Theory of Elasticity,” page 


230, equation 129, second edition, McGraw-Hill Book Company, Inc. 
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in the polar coordinate system. VY’ is called the Laplace or Har- 
monic Operator. Thus equation (2) may also be written: 


Assume a stress function such that y = y,+wy,,¥, will be the 
stress function valid for the case if there is no hole present. % » 
will “correct the error committed” by not considering the hole. 

The following boundary conditions have to be satisfied: 


M , 
» t. * eee ~~ at r=x-—- © (well known equation for 
an intact bent plate) 
b) Trr = 0 at r=a (no radial stress exists 
at the hole) 
c) T9790 at r=a (no tangential stress ex- 


ists at the hole) 


hence from Eq. (la): 
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y =rsin@ and sin’@ = +(3 sin@ - sin(36)], thus y becomes 
_ Mr ' , 
q& *- 94, (2 sin® - sin(3@)] _ «= 6 = © 2's «  @& 


It may be shown that the stress function - satisfies the bi- 
harmonic equation: vty, = 0, thus the stress system shown be- 
low is compatible. 


The stress components: 











1 @ 1 3? M . = 
i, “2 a aS = - 77 [siné + sin(3 4)] a, 
2 
96, = EE sing - sins 9) oP 
1 . 
ar ss « Cs ) = Tp [cose - cos (3 @)] c, 


It is convenient to assume ¥, in a similar manner as: (see 
Eq. 4) 
y 


+. f, (r)sin@ + f(r) sin (3 6) 


According to Goursat’s Theorem a Biharmonic Function (such 
as y, ) may always be expressed in the following form: 


oe r?G where YV?F =0 and V?G =0 


in other words F and G are harmonic functions. Since f (r) sin@ 
also has to be a biharmonic function it can be assumed that it is 
composed of a harmonic term and the product of another harmonic 
term and r’. If f,(r)sin@ were harmonic it would satisfy Laplace’s 
equation. 


Hence: 


02 1 


a 2 
Far [f (r) sine} + r Or Uh (Fr) sin] + 2 [f, (r) sing] =0 


ee 4 
r’? a0 
which may be rewritten as follows. 


siné 
r2 


, 3 sin 6 
sind - [f,(r)] + 








C) , 
sf, - SB" tf, @) = 0 
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simplifying by siné 


1 2 1 ' 
or? 4) + oop 8) oa 4,00,2.9 «2. 9 Me 2 3 @ 


The solution of Eq. 6 is the following: 
f (x) = A’r + Br= 


f, (r) is, however, biharmonic hence another term (Cr + Dr™)r’ 
has to be added. For A’ + D=A one obtains 


1 


(eyeAre Be 2+Cr -l& « © 8 S cit te of} 


The expression f,(r) may be found in a similar way. It can be 
shown that: 


f(r) = Er’ +Fr°+Gr*+Hr™ « ewele hf « « @ 
Thus y, =f, (r)sin@ + f, (r) sin(34)is now a biharmonic expres- 
sion and a stress system derived from it by means of equations 7, 
8 and equation 5 will be compatible and in equilibrium. 
The unknown constants: A, B, C, E, F, G and H may be deter- 
mined by satisfying the boundary conditions. 








y = (Ar + Br“ + Cr )siné 
(Er®+ Fr*+Gr™"+ Hr )sin(39) - - - - (9) 
Then 
’ lay, . 1 a*y, 
‘Tr, r or r* 3a6* 
os + 3Cr) sin + (sE° © err «3 « 3H) in (3 6) 
A,B - ; = = 
‘és (s ; Cr) sin - (Er Fr 7 73) 9sin (3 }) (10) 
at x ~ co T = fT 
xx rr 
But ee “ - Ter, thus Trr, =0 at x=r- co 


Therefore, the sum of the terms shown on the right hand side of 
Eq. 10 equal zero. 
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This is only possible if C = E = F = 0 








Hence Trr, 2% 2% sine - Che + o#) sin(3@) - - -(11) 

Constant A vanishes in "8, and To99 also since 

0,” : (3) FF Te + 
thus Eq. 9 becomes: 

vy, = Br™ sind + (Gr™ + Hr™)sin(3@) - - - - - -(12) 
From Eq. 12 (replacing y, by y, in Eq. 5) 

tT.» * 2Br~cosé + (2Gr™+ 4Hr™) 3cos (36) ae one 
and 

T =2Br~sing + (2Gr*%+12Hr™)sin(36@) - - - -(14) 


Pus 7 Re = 0 at r=a 


Using equation (5a) and 11 and substituting “a” for “r” one may 
find that: 





T = T + T = 0 at r=a 
It can be shown that this leads to the following values: 


_. Ma* _ Ma’ 
Gr-"o7 4. 8" i981 





Thus from Eq. (4) and (12) 


Mr’ ,.a° a‘ , 
41 (2-4 -3s1+ 1) sin (3 @) 


M 4 a 
wad +! es. + le + 
vyryuUtY, ri a“) sinég 
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Also: 


4 5 
Terr = Trr,* Trr,* mr (a, 1) sin @ + Ei 43+ Sa" - 1) sin (3 6) 


= Mr yy 3) coso+ MF (43, 32,..1) cos (3 4) = 
‘ré r 6, r0.— 4I ~,1/ cos r- r*” cos } 


a® a‘ 


_. Mr Mr , , 
T99* Too, * 99,7 - 1 ** 3) sin@ + 41 (4-3 - rit sin (34), 





Equations (15) determine the three stress components at any 
point of the plane given by r and ¢@. These stresses are compatible, 
while satisfying the equilibrium, and boundary conditions. The max- 
imum occurs at the circumference of the hole. (r = a) 


Here Trr th al. 


99) = . M4 [sing - sin(36)}- - - - - - -(16) 


It is known that: 


sin(3@) = 3 sin@ - 4 sin*é 
Therefore: 

: - Ma,,, 12 

(Tp) he seeeeare «Rs «© - os 
r=a 

(Tp) = 0 for 6 =n 
za 

or 6 * (2n + 1) where n = 0,1,2,3, --- 


In order to find the maximum value of (7 differentiate it 


@@ r=a 
with respect to @ and set it equal to zero. 
d 
ae = Ma i192 sin?@ cos@ - 2 cosé@) = 0 
dé I 
Thus 


12 sin?9@ = 2 Hence 6 =n7a * 24.15° 
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or if cosé@ = 0 so siné =+t1 and 
6 = ; (2n + 1) 


where again n = 0,1,2,3, --- 


The latter value yields a maximum stress of 


ny 


Kea 


Hence the maximum stress is twice the stress calculated by 
means of the elementary theory. In other words the factor of stress 





concentration is 2. 
It is interesting to note that (766)... reaches the maximum 


value at @ = ; (2n +1) but it is not necessarily the maximum stress 
of the whole system. 
Denoting 7 = . and 6 = ; (2n +1) in Eq. (15) one may obtain the 


following: 





Thus 





796 ~ I (18) 


In order to draw the stress diagram, let a = 1 and r vary from 
1 to 7. From Eq. (18) it is easy to establish a chart and the stress 
diagram shown below. (See Fig. 3) 





a=1 
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The dotted line represents the 


stress distribution without the pl 

hole, and the solid line indi- eG 

cates the present one. cc 

s( 

Thus the effect of the hole will be negligible at a few diameters fo 

distance away. fc 
The negative sign denotes the compressive stress and the posi- 

tive sign denotes the tensile stress. (a 
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A Graph for Determining CSP-5 Continuous 
Sampling Plans 


By John S. White 
Mathematics Group 
General Motors Research Laboratories 


Two new single level continuous sampling plans were recently 
proposed by Derman, Johns, and Lieberman (1). These plans, called 
CSP-4 and CSP-5, are essentially modifications of Dodge’s original 
continuous sampling plan CSP-1 (2). The new plans are best de- 
scribed by setting forth their procedures as compared with that 
for CSP-1. The procedures for CSP-1, CSP-4, and CSP-5 are as 
follows. 


(a) At the outset, inspect 100% of the units consecutively as pro- 
duced and continue such inspection until i units in succession 
are found clear of defects. 


(b) When i units in succession are found clear of defects, discon- 
tinue 100% inspection, and inspect only a fraction f = 1/k of the 
units, 


(b-1) selecting individual sample units one at a time from the 
flow of product, in such a manner as to assure an un- 
biased sample. 


(b-4,5) choosing the unit to be inspected at random from a seg- 
ment of size k. 


(c) If a sample unit is found defective, 


(c-1) revert immediately to 100% inspection of succeeding 
units and continue until again i units in succession are 
found clear of defects, as in paragraph (a). 


(c-4) revert immediately to 100% inspection, eliminating from 
the production process the remaining (k-1) units in the 
segment, and commencing 100% inspection with the next 
unit following the eliminated segment. Continue 100% 
inspection until again i units in succession are found 
clear of defects, as in paragraph (a). 
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(c -5) screen the remaining (k-1) items in the segment. Upon 
completion of this screening, commence 100% inspec- 
tion with the next unit produced. Continue 100% in- 
spection until again i units in succession, not including 
the (k-1) screened items, are found clear of defects, as 
in paragraph (a). 


(d) Correct or replace with good units all defective units found. 


Each of the plans CSP-1, CSP-4, and CSP-5 follows the above 
procedure except that each follows the subsections of paragraphs (b) 
and (c) which correspond to the plan number. 

If the production process is “in control” with a defect rate of 
p(q = 1 - p) then the respective AOQ functions are 


k-1)pq' 
AOQ(CSP-1) = yx 


(k-1)pqit? 
1 + (k-1)q?*? 


AOQ(CSP-4) 


=. _ (k-1)pq'** 
AOQ(CSP-5) = 1 eat 

From the form of the AOQ functions we see that the AOQL of 
a CSP-4 plan with a given f = 1/k and i is the same as that for a 
CSP-1 plan with the same f but with clearance number i+1. Thus 
specific CSP-4 plans may easily be obtained from available tables 
and graphs (2,3,4) for CSP-1. 

Although the AOQ function for CSP-5 differs from that of CSP-1 
by only a factor of q, there appears to be no simple analytic method 
of obtaining the CSP-5 AOQL from that of CSP-1. The AOQL for 
CSP-5 plans for the set of sampling fractions used in H-107 was 
computed and the results appear in Figure 1. For i = 50, the exact 
AOQL’s for certain f values are given in Table 1. 
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AOQL IN PERCENT 
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Introduction to Nonlinear Tchebycheff 
Approximation 


By John R. Rice 


1. THE PROBLEM. 


The linear abproximation problem is as follows: Given a func- 
tion f(x) and a set of n functions {g; (x)}, find the parameters (or 
coefficients) a of the approximating functions 


(1.1) L(A,x) 


which makes L(A,x) as “close” to f(x) as possible. There are a va- 
riety” of ways to measure how “close” L(A,x) is to f(x). In this 
article we consider the Tchebycheff (or uniform) measure. Thus 
L(A,x) is required to minimize 


(1.2) max iL(A, x) - f£(x)| 
X € [0,1] 


The specific interval [0,1] of approximation is used throughout this 
article, though the discussion is valid for any interval. 

The nonlinear approximation problem is identical with the linear 
approximation problem except that the approximating function F(A,x) 
no longer depends linearly on the parameters. The approximating 
function F(A,x) is assumed to depend continuously on the parame- 
ters. The problem is formally stated as 


*The most common way is to minimize 
J [L. A,x) - f(x) }? dx 


For linear approximation this least squares measure has the advan- 
tage that the parameters may be computed directly. However, there 
are no direct methods for nonlinear least squares approximations 
and there is some evidence that nonlinear Tchebycheff approxima- 
tions may be computed with less difficulty for certain approximating 


functions. 











JOHN R. RICE 
NONLINEAR TCHEBYCHEFF PROBLEM: Given f(x) continuous 
on [0,1| determine parameters A* so that 


max |F(A*,x) - f(x)| = max |F(A,x) - f(x)| 
x € [0,1] x € [0,1] 


for all parameters A. 


This article gives an informal account of some recent develop- 
ments in the theory associated with this problem. The method of 
presentation is to first consider very simple examples in some de- 
tail. An “intuitive feel” for the problem is thus obtained. The gen- 
eral results are then obtained from an abstraction of the essential 
elements of these examples. No formal proofs are given. 

There are three principal questions associated with this prob- 
lem. They are 


1. Does a best approximation exist? 
2. If so, is it unique? 
3. What properties are characteristic of a best approximation? 


The first two questions are, of course, of basic interest, but the 
answer to the third question is the most important from the practi- 
cal viewpoint. For it is these properties which distinguish the best 
approximations from other approximations, and the methods of 
computation are based on them. 


2. APPROXIMATION BY POLYNOMIALS. 


First let us consider the simplest possible case, namely when 
F(A,x) depends on one parameter a and 


F(A,x) = a. 


For any approximation a, to f(x), consider the error curve f(x) -a 
in Figure 1. The maximum error occurs at x,. It is easy to see 
that if we increase a slightly (say to the dotted line), we can de- 
crease max |f(x) - aj. Indeed we can increase a until the error at x, 
becomes as large as the error at x,. At this point we cannot im- 
prove the approximation by varying a. If we decrease the error at 
xX, we increase it at x, and vice versa. Thus in this simple case we 
can get a formula for the best Tchebycheff approximation, 


(2.1) a* = ; {max f(x) - min f(x) | 


xe[0,1] xe([0,1] 








at 
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Figure 1. Approximation by a Constant 


We now consider the next simplest problem, namely approxima- 


tion by 
(2.2) F(A,x) = a + bx 


We again examine our error curve in Figures 2a, 2b, which repre- 
sent the same situation. 











Figure 2. The Error Curves of a Trial Approximation to f(x). 


The largest error occurs at x,. Clearly we can adjust the constant 
a, so as to decrease the error at x, at the expense of increasing it 
at x,. However, we can actually do more. We can reduce the error 
at both x, and x, by adjusting both of the parameters a and b. Such 
an adjustment is shown by the dotted line. This process may be 
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Figure 3, In order to reduce the error at three 


points the approximation must cross a* + b* x twice. 


continued until the error at x, becomes as large as the maximum 
error. As long as the three errors are not the same size, the 
largest two may be reduced by varying a and b. Once the three 
errors are equal we cannot improve the approximation. This is 
illustrated in Figure 3. Since the errors are all equal a better ap- 
proximation than a* + b* x must cross a* + b* x twice. This im- 
plies that the better approximation is curved, which is not possible 
for a line. Note that it is not sufficient for the maximum error to be 
assumed three times, but the maximum error must be assumed 
three times with alternating sign. This is illustrated in Figure 4. 





Figure 4, The maximum error is assumed 


4 times, but a + b x is not a best approximation. 
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In Figure 4 if one decreases the error at x, and x,, it is automati- 
cally decreased at x, and x,. 

We have come to the conclusion that a* + b* x can be a best ap- 
proximation only if the error curve appears as in Figure 3. We 
give this situation a name in the following definition: 


DEFINITION 1. The error curve F(A,x) — f(x) is said to alter- 
nate n times if these are n +1 points 


OS2,< 24,4 ...SBans 


1 2 


such thai 
(2.3) F(A,x;) - f(x; ) = -[F(A,x;,,) - £(x;,,)] = = max | F(a,x) - f(x)| 


We may restate this conclusion by saying: 


A necessary condition for a*+ b* x to be a best approximation 
to f(x) is that a* + b* x - f(x) alternates at least twice. 


We may also conclude from Figure 3 that 


A sufficient condition for a* + b* x to be a best approximation to 
f(x) is that a* + b* x - f(x) alternates at least twice. 


In these two statements we have answered the third question as- 
sociated with the approximation problem for (2.2). Best approxi- 
mations to f(x) are characterized as exactly those approximations 
a* + b* x such that a* + b* x - f(x) alternates at least twice. 

There is a slightly different way to look at approximation by 
a+bx. There are two parameters which may be varied in a + bx. 
They may be used to “push down” the error curve at any two points. 
When one “pushs down” the two largest “humps” on the error curve, 
there is a corresponding rise in other humps. As soon as a third 
hump of the error curve becomes as large as the two largest ones 
(including alternation of sign) the best approximation has been 
reached. The fact that the error curve may be “adjusted” at any two 
points is due to the fact that the parameters are independent, i.e. 


(2.4) 








for x, # x.. 
1 2 
We may briefly examine this line of reasoning for approximation 


by 
(2.5) F(A,x) = a + bx + cx’ 
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a 


Figure 5. Alternation 2 and 3 times for a quadratic approximation, 


The situation of alternation 3 times is illustrated in Figure 5. In 
Figure 5b a + bx + cx” - f(x) alternates only twice and we can find 
a parabola (i.e., a function of the form (2.5)) which decreases the 
error at each of the three points where the maximum error is at- 
tained. Such a parabola is shown by the dotted line. In order to re- 
duce the maximum error in Figure 5a an improved approximation 
(shown dotted) must intersect a* + b*x + c*x’ at least 3 times. This 
is impossible for parabolas. We have then arrived at the same type 
of characterization for quadratic approximation as for approxima- 
tion by a + bx. 

We now turn to the most general linear approximating function 
(1.1). There were two properties of polynomials that were used to 
obtain the characterization of best approximation by alternation. 
The first was the ability to “push down” the error curve at any n 
points (for n parameters). For the approximating function L(A,x) 
we may obtain this ability by assuming 


Q (&,) 9, &,) .-.. Q (x) 
(2.5) PY, (X,) ,(x,) ... Gy X,) #0 
Pon (K,) Pnlk.) .-- G(X) 


for any set of n distinct points {x; }. This condition then implies that 
alternation at least n times is necessary for best approximation. 
The second property of polynomials used is the fact that two distinct 
polynomials of nth degree can intersect at most n times.* This 








*This is more commonly stated by saying that a polynomial of 


nth degree has at most r zeros, 
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property for L(A,x) also follows from condition (2.5). Thus if we 
assume the set of functions {@ (x)} satisfy (2.5) we find that best ap- 
proximations are characterized by alternation of the error curve at 
ieast n times. A set of functions satisfying (2.5) is called a Tche- 
bycheff set [1], [2]. 

We now turn to the first two questions. It is always true that 
there is a sequence {A; |i = 1, 2, ...} of parameters so that 


Lim max |L(A; ,x) - f£(x)| = inf max |L(A,;,x) - f(x)| 
A 


li- © 


We need only to know that this sequence possesses a limit. Due to 
the linear dependence of L(A,x) on A, the set {L(A; .x)} is an equi- 
continuous set of functions. Thus the sequence {L(A;,x)|i = 1, 2,...} 
contains a convergent subsequence and the limit of this subsequence 
exists and is a best approximation L(A*,x) to f(x). 

The uniqueness of the best approximation by L(A,x) follows from 
the fact that L(A*, ,x) - L(A*,,x) can have at most n - 1 zeros. This 
is a consequence of (2.5). 

We can now state the following theorem which gives the princi- 
pal points of linear Tchebycheff approximation. [1], [2], [11]. 


THEOREM 1 - Let f(x) and 9; (x) be continuous on [0,1]. If 
{o;()} satisfies (2.5) then we have 


A. a best approximation L(A*,x) to f(x) exists 
B. L(A*,x) is unique 


C. L(A*,x) is the best approximation to f(x) if and only if 
L(A*,x) - f(x) alternates at least n times on |0,1) 


Theorem 1(A) follows from the fact that L(A,x) is linear in A. 
Theorem 1(B) follows from the fact that L(A, J) « L(A, ,x) may have 
at most n-1 zeros; the if portion of Theorem 1(C) also follows 
from this. The only if portion of Theorem 1(C) follows from the fact 
that we may “push down” the error curve at any n points. Actually, 
this theorem follows from only two things, linearity and (2.5), how- 
ever when we consider nonlinear approximation condition (2.5) will 
disappear and be replaced by conditions similar to those given here 
which are (in the linear case) consequences of (2.5).* 


~ *Motzkin [3] and Tornheim [10] have defined unisolvent families 
and n-parameter families, respectively where properties equivalent 
to the “n-l zero properties” and “push down” property are ab- 
stractly assumed. For any nonlinear approximating function with 
these properties Theorem 1ABC is valid. The elementary nonlinear 
approximating functions do not have these properties. 
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3. APPROXIMATION BY AB* + C ON [0,1]. 


We now consider a simple example of a nonlinear approximating 
function in order to see what one may expect. We consider approxi- 
mation as follows: 


(3.1) ab* + c ~ f(x) x € [0,1] 


The first question that we need to consider is: do best approxima- 
tions exist for all continuous functions? If we have a sequence 
{a;b;* + ¢;} so that 


(3.2) Lim max la; b; *+c; - f(x)|= inf max lab*+c - f(x)|, 
i oo ’ a,b,c 


this sequence must certainly remain bounded. Some possibilities 
for behavior of this sequence are shown in Figure 6. 











(b) (c) 


Figure 6. Some exponentials. 


The limits for case (a) may be 


a+c x =0 
(3.3) Dix) ={ 
c 0<x31 


However, this approximation is no better than the constant c. This 
is a manifestation of the fact that one cannot improve an approxima- 
tion by changing its value at only one point. Hence we identify all 
“curves” (3.3) with the constant c. The same approach applies for 
Figure 6c. In Figure 6b, nothing much happens; the exponentials just 
settle down to another exponential OR to a line. A little thought 
shows that we cannot find a best approximation to a straight line. 
The maximum error may be made arbitrarily small, but a b*+c is 
never actually a line. We may avoid this embarrassing situation by 
simply extending the definition of our approximating functions to 


ab*+ec b #1 
(3.4) E(A,x) -{ 


ec +dx b=lora=0 
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For L(A,x) we are able to “push down” the error curve at any n 
points (here to number of parameters is 3, so n=3). Let us ex- 





” amine E(A,x) for a similar property. See Figure 7. 
L- 
e 
'S 
| 7— "2 
: Figure 7. Improvement of the approximation at 3 points. 


To approve the approximation at x,, x, and x, one increases c 
slightly, decreases a and increases b. One can do the same thing 
for a straight line. Thus it is seen that if the error curve alternates 
only twice then one can improve the approximation. However, there 
is one exception and that occurs where b = 1, d = 0 (or alternatively 
a=0). 





eee 


Figure 8. The degenerate case of a b* + c = constant. 


There is no possibility for a function E(A,x) to behave as the 
dotted line in Figure 8, for every E(A,x) is monotonic. Hence, we 
arrive at the following conclusion: 
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If b# 1, d# 0 (or a#0O) then 3 alternations of the error curve 
are necessary for E(A*,x) to be a best abproximation to f(x). If 
b = 1, d = 0 (or a = O) then only 2 alternations are necessary. 


Thus the number of alternatives necessary for a heat approximation 
depends on the best approximation. The question of interest at this 
point is: does the same result hold for the number of alternations 
sufficient for a best approximation? Recall from the linear theory 
that this question relates to the number of zeros that the difference 
of two approximating functions can have. It is easy to show that 


(3.5) E(A, ,x) - E(A, x) 


can have at most 2 zeros, and if E(A, ,x) = constant then (3.5) can 
only have one zero. Therefore we have the following. 


CHARACTERIZATION. E(A*,x) is a best approximation to f(x) 
if and only if E(A*,x) - f(x) alternates k times where k = 2 if 
b= 1,d= 0 (ora =0) and k = 3 otherwise. 


The uniqueness of the best approximation follows from the dis- 
cussion of the number of zeros of (3.5) in the same way as for linear 
approximating functions. 

It should be noted here that a condition analogous to (2.5) is not 
present. The questions of the number of zeros of (3.5) and the abil- 
ity to “push down” the error curve have become separated. The 
question of the existence of best approximations has remained a 
separate problem. 


4. ABSTRACTION TO VARISOLVENT FUNCTIONS. 


In this section we extend the theory indicated by the example of 
a b* + c to an abstract approximating function F(A,x). It is assumed 
that F(A,x) depends continuously on A and x. The domain of param- 
eters A is denoted by P. The set P is viewed only as a set to index 
the various functions represented by F(A,x). It is sometimes possi- 
ble (but often not easy) to identify P with a simple subset of Euclid- 
ean space of appropriate dimension. 

In the preceding section we have seen that there is a character- 
istic number of alternations for best approximation. This number 
is not a constant, as for linear approximating functions, but rather 
is a function of the particular best approximation, or equivalently, 
of the parameters A. This idea can be abstracted as follows: 


DEFINITION 2. F/(A,x) is a varisolvent function if there exists a 
degree m(A) defined on P such that F(A*,x) is a best abproxima- 
tion to f(x) if and only if F(A*,x) - f(x) alternates m(A*) times. 
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This is not the usual or original definition of a varisolvent func - 
tion. The original definition [6] says F is varisolvent if F has the 
following properties: 


PROPERTY Z. F(A,x) - F(A*,x) has at most m(A*) - 1 zeros. 


PROPERTY OF SOLVENCE. given {x;|j = 1, 2,..., m(A*)} 
and «> 0 then there existsa 5(A*, €, x 4 > 0 such that 


|F(A*, x;) - Y; I< € 
implies that there exists a solution A to the equations 
F(A, x;) = Y; i «3, @ <.., S94 
with max \F(A,x) - F(A*,x)|< e. 


Property Z is the property of linear approximating functions 
used to establish the sufficiency of m(A*) alternations for best ap- 
proximation. This property also implies that best approximations 
are unique. The Property of Solvence is an « - 5 statement of what 
is meant by being able to “push down” the error curve at any m(A*) 
points. It is one of the deeper results [5] of approximation theory 
that these two definitions are equivalent.* 

For varisolvent approximating functions the characterization 
and uniqueness of best approximations are intimately connected by 
Property Z. The question of the existence of best approximations is 
completely separated. 


5. EXAMPLES. 


In this section a few known and conjectured examples of vari- 
solvent functions are given. 


1. Rational Function {1}. 
Dax , s(x) > 0, x € [0,1] 


(5.1) R(A,x) = s(x) 
mb -x/ 


bx > 0, x € [0,1] 


, }-0 


Ua? +b? =1 


1 


*In order to show equivalence one must assume that P is anarc- 
wise connected set and that {F(A,x)| max |F(A,x)| = M Sis closed 
under pointwise convergence, 
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This is the classical example of a varisolvent function. Let 

(5.2) P, ={Ala®=a® =... sank = pM = pM =... = pm-K =} 
If Ae Py and A¢ P,.,, then the degree of A is 

(5.3) m(A) =n+m-k +1 


Normally we have k = 0. 
2. Exponential functions {7}. 


7 
(5.4) E(A,x) = D | 


my) 
i=1 J 


D aij x! | ei lagi, 1 ty 1< © 


k 
(5.5) 2 (mj + 1) =n 


1=1 


This is the approximating function one obtains when the defini- 
< tix 
tion of the function »D a;e* is extended to include the polynomial 
1=1 


limits. The degree of Ais 


x 

(5.6) m(A) =n + 2 (m; + 1) 
1=1 

3. Constrained linear functions [6}. 
n 

(5.7) L.(A,x) = 2 a, 9; (x) 


i=1 


where the ith parameter a; is constrained to lie in the open interval 
Ij. Condition (2.5) is assumed. Then the degree of A is 


(5.8) m(A) =n 


4. Exponential Cosine (Conjectured). 


(5.9) E.(A,x) = ae™ cos (@ +x) +b 





(5.1 


(5.1 
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5. Logarithmic (Conjectured), x ¢€ [-1,1]. 


n 


(5.10) Log (A,x) = 2) a; log (1+b,x) 


_e | . b; 6 ¢ 1 

6. Exponential ratio (Conjectured). 

n 

> a; etix 

1=1 
(5.11) eh ieee la;l,i bs] < ec 

m 

27 b; e5i* It; l,1s;1 < © 

1=1 


6. COMPUTATION METHODS. 


At the present time only a few computational methods are known 
to apply to these examples. See [4], [9] for some discussion of 
these methods. There is a conjecture however that the Remez algo- 
rithm [9] may be modified so as to apply to all varisolvent func- 
tions. The basis idea of this method is illustrated in Figure 9. 








Figure 9. The Remez Algorithm for m(A,) = 3. 


One is given an estimate F(A,,x) of the best approximation to f(x). 
It is required that the error curve F(A;,x) - f(x) alternates m(A;) 
times. If it does not one chooses the m(A;,) largest local extrema 
of the error curve which have alternating signs. These occur at x , 
x, and x, in Figure 9. One then finds the estimate F(A,4,,x) to the 


best approximation F(A*,x) which alternates on these points i.e., X.» 
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x, and x,. This approximation is shown in dotted line. It is known 
that this algorithm converges for linear approximating functions and 
some nonlinear ones [4], [8], [9]. It is conjectured that it converges 
for all varisolvent approximating functions. 

There is one essential modification that must be made to this 
algorithm due to the varisolvence of the approximating functions. 
If F(A;,x) is of degree, say 5, it is possible that the degree of 
F(A;,,,x) may be less, say 3. This implies that one cannot find any 
F(A,x) that alternates 5 times on the six local extrema of F(A; ,x) 
- {(x). This will, of course, become apparent during the process 
used to compute F(A;,,,X). In such a case, one knows that the de- 
gree of F(A,;,,,x) is less than 5 and hence one replaces the 6 ex- 
trema by 5 extrema with largest size and with alternation. One then 
again attempts to compute F(A;,,,x). One proceeds in reducing the 
number of extrema used until F(A,;,,,x) is found. 
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